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00 ■ Abstract. We present preliminary results of the application of a new sophisticated 

£Nj ■ code which allows high precision integration of orbits of stars belonging to a dense 

stellar system moving in the vicinity of a massive black hole. This mimics the situation 
observed in the center of many galaxies, where a nuclear star cluster contains a massive 
■ black hole which, in the past, was, likely, an active engine of violent emission of radia- 

| tion. The main scope of our work is the investigation of the relaxation of the super star 

cluster on a sufficiently long time, together with the investigation of its feedback with 
' the massive black hole. 
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1. Introduction 

It is known, nowadays, that many galaxies harbor a supermassive black hole (SBH) 
in their innermost regions. In many cases, the massive black hole is surrounded by 
an ultra compact stellar system, the so called Nuclear Star Cluster (NSC). Only in the 
case of the Milky Way the NSC is resolved in its stellar components, thanks to the high 
. sensitivity and spatial resolution of modern instrume nts. The formation of c ompact 

nuclei in galaxies has been studied, among others, by ICapuzzo-Dolcettal (1993) which 
developed the original iTremaine et al.l(ll975l) hypothesis of the M 31 nucleus formation 
via repeated infalls of globular clusters braked by dynamical friction, as alternative to 
the in-situ form ation hypothesis. Nevertheless , the fu ll jV-body simulations performed 
so far (see, e.g.. lCapuzzo-Dolcetta & Miocchil d2008h . lAntonini et all d201ll) ) have not 



been extended enough in time and with a sufficient resolution to reach fairly conclusive 
statements about it. 



2. Our new N-body code 

To the purposes discussed above, we developed a direct summation Af-body code which 
includes a Hermite 6th order time integrator implemented by mean of a hierarchi- 
cal block time-stepping. This code is written in C++, and uses MPI, OpenMP and 
CUDA to fully exploit the available computational power of modern hybrid compu- 
tational clusters. Our first use of this new Af-body code is to studying the dynamical 
evolution of the Milky Way NSC over few relaxation times. 

In modeling the NSC, we take as initial conditions a mass Mjvsc - 1-32 x 10 7 M Q 
distributed over N 270, 000 equal mass bodies, and a central SBH of mass Msbh - 
4 x 10 6 M Q . All the preliminary tests have been executed on our hybrid machine com- 
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posed by 2 NVIDIA Tesla C2050 GPUs in a host with 2 Intel Xeon X5650 esacore 
CPUs. As we see in Fig. [T] the improvement of performance is crucial when the 
number of particles to be updated is small because, in this case, the time spent in trans- 
ferring data among CPUs and GPUs becomes comparable to that required by the force 
calculation (the GPU kernel). Moreover, we verified that the update frequency is much 
higher for small numbers of bodies (as expected given the presence of an SBH) and this 
constitutes the real bottleneck for such kind of simulations. 
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Figure 1 . Time spent in the execution of different parts of our code, in function of 
the number of particles to be time-updated. 



3. Final remarks 

The simulation of the long term evolution of a Super Star Cluster around a massive 
black hole is a real computational challenge. At the light of the analysis presented here, 
we believe that a hybrid cluster of intermediate size fully exploiting the computational 
efficiency of our Hermite's 6th order code will allow to reach the goal of the description 
of the long term evolution of dense super star cluster around a massive black hole. 
Actually, we estimate that a sufficiently extended (in time) simulation will require about 
100 days of computation per node, as extrapolated from the preliminary benchmark 
runs done on our private machine. 



References 

Antonini, E, Capuzzo-Dolcetta, R., Mastrobuono-Battisti, A., & Merritt, D. 2011. Submitted, 

see 2011arXivlll0.5937A 
Capuzzo-Dolcetta, R. 1993, ApJ, 415, 616 
Capuzzo-Dolcetta, R., & Miocchi, P. 2008, MNRAS, 388, L69 
Tremaine, S. D., Ostriker, J. P., & Spitzer, L., Jr. 1975, ApJ, 196, 407 



